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O ■ Abstract 

■ An algorithm which computes a solution of a set optimization problem 
1 is provided. The graph of the objective map is assumed to be given by 

finitely many linear inequalities. A solution is understood to be a set of 
points in the domain satisfying two conditions: the attainment of the in- 
' fimum and minimality with respect to a set relation. In the first phase of 

, the algorithm, a linear vector optimization problem, called the vectorial 

relaxation, is solved. The resulting pre-solution yields the attainment of 
the infimum but, in general, not minimality. In the second phase of the 

■ algorithm, minimality is established by solving certain linear programs in 

combination with vertex enumeration of some values of the objective map. 
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OV 

^ 1 Introduction 

In this work we consider the problem to minimize a set- valued map F 

(N 



with polyhedral convex graph with respect to the relation 

F(x)<F(u) : F(x) + C 2 F(u), 



where C denotes a polyhedral convex ordering cone that contains no lines and 
has nonempty interior. The objective map can be considered as a function from 
R" into the space Q of all closed convex subsets of M. q . With the above ordering 
relation, one obtains a complete lattice, i.e., the infimum inf {F(x) \ x e W 1 } in 
the sense of a greatest lower bound with respect to < always exists. Solution 
concepts for complete- lattice- valued problems have been introduced in [T7j . The 
main idea is that, beyond scalar optimization, minimality and infimum attain- 
ment arc two different conditions and a solution shall involve both. Such a 
solution concept is also useful in vector optimization J55J H3 [201 HI] ■ In a set- 
valued framework, it has been used, for instance, in [51 [11] [T5] . Applications of 
set optimization based on the above ordering relation and solution concept can 
be found in mathematical finance in the framework of markets with frictions, 
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see e.g. [T31 HI QUI Q31 EH H2J- But, in specific calculations, only the infimum 
attainment have been considered so far. The algorithm presented in this note 
ensures also minimality. 

Optimization problems with set- valued objective function but partially based 
on other ordering relations or other solution concepts have been investigated by 
many authors. References and results can be found, for instance, in [331 El 113 

ei s cei ess eng. 



2 Preliminaries 

A set P £ W is said to be polyhedral convex if there is a representation 

r 

P=0{yeRi\(w l ) T y>^} (1) 
i=i 

where w 1 ,...,w r e R 9 \ {0} and 71,..., 7 r e R. Equation (QJ is called H- 
representation of P. Every non-empty polyhedral convex set PcRi can be ex- 
pressed as a (generalized) convex hull of finitely many points x 1 , . . . , x s e R 9 (s e 
{1,2,3,...}) and finitely many directions d\...,d* e R«\{0} (t e {0,1,2,...}) 
through 

P = I J] AiX* + 2 M j A, 5= 0, J] A( = 1, Hj > 1 , (2) 

U=l 3 = 1 i=l j 

where d e R 9 \{0} is called a direction of P if P + {A • d} c P for all A > 
0. This can be also written with the convex hull of the points and the cone 
generated by the directions as P = conv {x 1 , . . . , x 3 } + cone {d 1 , . . . , d'}. We 
set cone0 = {0}, thus, P is bounded if and only if t = 0. Equation Q is called 
a V -representation of P. Numerical methods to compute a V-representation 
from an H-representation and vise versa are called vertex enumeration, see e.g. 
[HE]- We denote by clP and intP, respectively, the closure and interior of a 
set P c W. 

We assume throughout that C R 9 is a pointed (i.e., C n (— C) = {0}) 
polyhedral convex cone with int C # 0. The cone C yields a partial or- 
dering on R 9 where y <c w is defined by v — y e C. If C = Ri := 
{y e R 9 | j/i > 0, . . . , y q 5= 0}, the component-wise ordering is abbreviated 

to <. The polar cone of C is the set C° := {v £ R 9 | Vy e C : v T y < 0}. A point 
y e R« is said to be C -minimal in P E R« if y e P and ({y} - C\{0}) n P = 0. 
We assume that an H-representation of C is given, that is, a matrix Z e R 9Xp 
such that 

C = { 2 /eR'|Z T y>0}. (3) 

Let 5 denote the family of all closed convex subsets of R 9 . By Qc we denote the 
subfamily of those elements P of Q having the additional property P = P + C. 
For P, Q e Q we define 

P<Q : P + C^Q + C, 

which can be equivalcntly written as P + C 2 Q. The ordering < is reflexive and 
transitive in Q (quasi ordering) and, additionally, antisymmetric in Qc (partial 
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ordering) . For P,Q e Q we define an equivalence relation by 
p ~q :<= ^. P + C = Q + C. 

Clearly, the quotient space Q/ ~ is isomorphic to Gc an d thus < is a partial 
ordering in Q/ ~. 

Remark 2.1. In a theoretical framework the space Qc is often more convenient 
and leads to easier formulations. From a computational viewpoint, however, the 
usage of Q and Q ' / ' ~ seems to be more natural. This is due to the fact, that an 
Pi /V -representation of some P £ Q with P ci P + C might be known whereas 
getting an H/V -representation of P + C would require computational effort. 

The partially ordered set (G/~, <) provides a complete lattice, i.e., for every 
subset of Q/~ there exist the infimum and supremum, see e.g. |20j for more 
details. To simplify the notation we express the infimum and supremum in terms 
of the (quasi-ordered) space <), where we have in mind that we actually deal 
with representatives of equivalence classes. Thus, for nonempty sets V c Q we 
have 

inf V = clconv (J (P + C) supP = f] (P + C). 

PeV PeV 

Furthermore, we set inf = and sup0 = W. To express minimality we 
define for P, Q e Q: 

P SQ ■ (P < Q and P 7^ Q). 

Let F : R" zz£ M. q be a polyhedral convex set-valued map, that is, its graph 

grP := {(x,y) e M" x R q \ y e F(x)} 

is a polyhedral convex set. We assume throughout that an H-representation of 
grP is known. This means, there are A e W nxn , B e R mxq and b e K m such 
that 

grP := {(x,y) el" x R q \ Ax + By > b} . (4) 

The domain of P is the set domP := {x e K n | P(x) # 0}. We consider the 
following set optimization problem: 

minimize P : E" =t R 9 with respect to < . (P) 

Problem <JPj) is called feasible if domP 0. We assume throughout that JFJ) 
is bounded in the sense that 

3w e W : {v} < inf Fix). 

aneR" 

In our setting, it is sufficient for (|P]) being bounded that domP is a bounded 
set and 

Vx6R ,l ,3v6R 9 : {v} < F(x), 

where the latter condition is obviously satisfied for a map P with bounded values 
F(x). Note that the algorithm introduced below can verify whether the problem 
is bounded or not. 

The following solution concept is based on a combination of minimality and 
infimum attainment as these notions do no longer coincide in vector and set 
optimization. It is an adaptation of the concepts introduced in jTTJ [501 H2] to 
the present setting. 



3 



Definition 2.2. A point x e domi* 1 is said to be a minimizer for ([P]) if there is 
no x e R™ with F(x) S F(x). A finite set X c dovnF is called a finite infimizer 
for (|Pj if the infimum is attained in X , that is, 

inf Fix) = inf F(x). 

A finite infimizer X of (jP|) is called a solution to <JPj if it consists of only 
minimizer s. 



3 Vectorial relaxation and pre-solution 



Consider the linear function / : R™ xR'^ M. q , f(x, y) = y. Because of formal 
reasons we understand / as a set- valued map whose values are singleton sets, 
i.e., 

/:»"xR'-tR ? , f(x,y) = {y}. 
The vectorial relaxation of the set optimization problem ([P]) is defined as: 

minimize / : 1" x 1' 4 R g with respect to < subject to y e F(x). (VR) 

Of course, (|VR[) can be seen as a special case of a set optimization problem, 
whence the above definitions apply also to (|VR|) . Obviously, (|VR[) is feasible if 
and only if so is ([P]). As / is single-valued and the constraint y e F(x) can be 
expressed by finitely many linear inequalities, (|VR[) is (equivalent to) a linear 
vector optimization problem. We have 



inf F{x) 



inf inf \y\ 

ssK" yeF(x) 



inf f(x,y), 

aeR™ yeF(x) 



(5) 



i.e., fP]) and (|VR[) have the same infima. This implies that |F]) is bounded if 
and only if so is (jVRI) . Equation ([S]) motivates the following concept. 



Definition 3.1. A finite set 
of © if there exist y i e W, 



G M 

1,.. 



, k} is called a pre-solution 



,,k such that {{x\y l ) T e R n x E 9 | i = 
1, . . . , k} is a solution of the vectorial relaxation (|VRj) of (|P|) . 

Of course, every solution is also a pre-solution. The following example shows 
that the x-component of a minimizer of (|VR|) is in general not a minimizer of 
()Pl) . This means that a pre-solution of (JPj is in general not a solution of ((Pj) . 



Example 3.2. Consider the set-valued map F : 
flU), gr F c R 4 is given by 



where, according to 
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Setting 



x 1 := (-1,0) T , x 2 := (-1,1) T , x 3 := (0,0), x 4 := (0, 1) T , 
we see that 

Fix 1 ) =conv {(0,2) T ,(2,0) T ,(0,3) T ,(3,0) T } ; 
^(* 2 )=conv{(0,2) T ,(-l,4) T ,(l,2) T }, 
F(x 3 ) =conv {(0,2) T ,(2,0) T ,(0,3) T ,(4,-1) T }, 
F(x 4 ) = conv {(0, 2) T , (2, 0) T , (-1, 4) T , (3, 0) T j . 

Consider the problems fF]) and (|VRp for the ordering cone C := R+. Set 

y 1 ■= (0,2) r , y 2 := (2, Of, y 3 := (-1,4), y 4 := (4, -if . 
The infimum for both problems ([P]) and (jVRp can be expressed as 

inf F(x) = conv {y 1 , y 2 , y 3 , y 4 } + R 2 + =: Q, 

set 2 

where each of the points y 1 ,...,^ 4 is ^c-minimal in the set Q. It follows 
that, for instance, the set {(x 1 ,y 1 ) T , (x 1 ,y 2 ) T , (x 2 ,y 3 ) T , (x 3 ,y 4 ) T } is a solution 
of (|VR|) . Hence {x 4 ,a; 2 ,a; 3 } is a pre-solution of (|P]). But {x 1 , cc 2 , x 3 } is not 
a solution of (|P]). Indeed, we have F(x 3 ) S F(x x ) and i^(a; 4 ) S F(x 2 ) which 
means x 1 and x 2 are not minimizers for problem (jPj). Note further that {x 3 , x 4 } 
is a solution to (|Pl) . 



4 Algorithm 

According to dU and d3]), let an instance of problem fF]) be given by A e K mxra , 

5 e K mXf ', b e R"\ Zel^. The algorithm computes a solution to (JFJ) if the 
problem is feasible and bounded. Otherwise it detects whether ([P]) is infeasible 
or unbounded. 

In the first phase of the algorithm, the vectorial relaxation (|VRj) , which is 
(equivalent to) a linear vector optimization problem, is solved. A solution can 
be obtained, for instance, with Benson's algorithm, see e.g. [21 [51 12"4"1 [231 [2U1 [T2] . 
We know that (|P]) is bounded if and only if so is (jVRp . But Benson's algorithm 
is able to detect if (jVRp is unbounded. Note further that in [20], Benson's 
algorithm was extended for unbounded linear vector optimization problems. 

In the second phase, for every point x° of the pre-solution obtained in the 
first phase, we construct a sequence (a; , x 1 , x 2 , . . . ,x ) with F(x°) ^ Fix 1 ) :> 
F(x 2 ) :> ■ • • F{x l ) until, after finitely many steps, a minimizcr x l is obtained. 

For parameters w e K ? and S e domF, we consider the following scalar 
problem: 

maximize w T y subject to y e F(x), F(x) < F(x). (P(w,x)) 

As the yi, . . . , y q e R are considered to be auxiliary variables, we use the follow- 
ing convention: x is said to be a solution to (P(w,x) ) if there exists y e M. g such 
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that (x, y) is a solution to (P{w,x) ) in the ordinary sense. In practice this means 
that (x, y) is computed but only x is used. Obviously, we have the following 
lower bound f3 for the optimal value a of (P{w,x) ): 

a(w,x) := sup {w T y\ y e F(x), F(x) < F(x)} 
^ sup {w T y\ y e F(x)\ =: j3(w, x). 

This leads to an optimality condition. 

Lemma 4.1. Let (|Pj be feasible and bounded. For some x e domf 1 let an 
H-representation of the set F(x) + C be given, that is, 

F(x) + C = {y e R«| (t^) T |, < lh j = 1, . . . ,r) . 

If arw 3 ,x) = P{w^ ,x) for all j e {1, . . . , r}, i/ien x is a minimizer for (jPj. 

Proof. Assume that a; is not a minimizer for ([P|), i.e., there exists x e R™ with 
^ f (S). Hence there is some y e -F(a:) and some c e C such that y + c ^ 
F(x) + C. Since C + C = C, we conclude that y $ F(x) + C . Thus there is some 
j e {1, . . . , r} such that w 3 y > "fj which implies aOw^ , x) > jj /3(tf J , x) ■ O 

With the aid of a V-representation of the set -F(ai) + C, that is, 

F(x)+C = conv {y\...,y s } + C, (7) 



(P(w,x)) can be transformed into a linear program. Note that the assumption 
that (|P[) is bounded was used in {?]), otherwise the cone on the right hand side 
can be a superset of C. By Q, the first constraint y e F(x) can be expressed 
as Ax + By > 6. The second constraint can be transformed as follows: 



F(x) < F(x) F(x) + C 2 F(x) + C 

Mi = l,...,s: i/'ef(i) + C 
<^> Vi = l,...,s, 3c l e C : ?/ 1 -c ?; e F(x) 

Vi = l,...,s, 3c 4 e C : At - Be 1 > b - By 1 



(8) 



Thus, (P(u>,x) I is equivalent to the linear program 
f Ax + By 5* b 

maxw T ys.tJ Ax - Be 1 ^ b - By 1 (i = 1, . . . , s) ^ (P(iu; y 1 , . . . , y s )) 
[ Z T c l > (i = l,...,s) 

which has n + q(s + 1) variables and m + ms +ps constraints. According to the 
above convention we will speak about a solution x e R™ and do not mention the 
q(s + 1) auxiliary variables y, c 1 , . . . , c s e R 9 . 



Algorithm SetOpt. 

Input: 

H-representation of gr.F according to Q: A e R mxn , B e R mx «, 6 e R 
H-representation of the ordering cone according to ([21): Z e R ?xp ; 
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Output: 

A solution X of jF]) if (JPj) is feasible and bounded, X = otherwise; 
The solution status for (|P|); 

Phase 1: 

solve (EE}; 

if (|VR[) is infcasiblc then status <— "jP|) is infeasible."; stop; end; 
if (IVRj) is unbounded then status <— "([0 is unbounded."; stop; end; 
store a pre-solution {a; 1 , . . . , cc fc } of (fPjl ; 
Phase 2: 

for i <— 1 to fc do 
flag - 1; 

while flag = 1 do 

compute a V-representation and an H- representation of F{x l ) + C: 
F(x l ) + C =conv {y 1 ,...,?/ s } + C 

= {yeW>\ (wtfy^v, j = l,...,r}; 
for j <— 1 to r do 

if u^/lw;- 7 1| ^ K then 
K^Kvj {w7|M|}; 

solve (P(w : ' ; y , . . . 2/ s )); £j <— solution; a <— optimal value; 
(3 <- max|(^)V, ■ ■ ■ , (w J ) T y s j; 
if a > /3 then break (i.e., exit the inner-most loop); 
end; 

if j = r then flag <— 0; 
end; 
end; 

X <- X u {^}; 
end; 

status <— "((Pj) /ias oeen solved."; 

We next show that the algorithm works correctly and is finite. We prepare the 
theorem by two lemmas. 

Lemma 4.2. Let (|Pj fee feasible and bounded. Consider some x e domf , a 
halfspace H = {i/el'l w T y < 7} containing the set F(x) + C and finitely many 
points y , . . . , y s e M? such that ((JJ) holds. Then, 

(i) The linear program (P(u>; y , . . . n s )) /ias an optimal solution; 

(ii) The lower bound (3 defined in © can oe expressed as 

(3(w, x) = max {w T y 1 , . . . , w T y s } . 

Proof. Let y e F(x). Then, Ax + By ^ b. By J7J), for all i e {1, . . . ,s}, we 
have y % e + C, i.e., there is some c 1 e C such that y l — c % e F(x), or 

equ ivalently, Ax — B e 1 5= b — By 1 . Hence, the point (x, y, c 1 , . . . , c s ) is feasible 
for (P(v);y\...,y'% 
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Since (0 is assumed to be bounded, there exists wel' such that {v} + C 2 
F(x) for all x £ R n . As H contains F(x) + C, we must have w £ C°. It follows 
that 

sup {w T y\ x e R", y e F(x)} w T v + sup {w; T c| c e C} = w T v, 



which implies that (P(u>,x)| and hence (P(u>; y , . . . , y s ) I is bounded. This 
proves (i). 

Statement (ii) follows from ((7]) taking into account that sup ceC w T c = for 
w e C°. □ 

Lemma 4.3. Let (jP| &e feasible and bounded and let x £ dom_F. If x e R™ is a 
solution of (P(w,x) ) for some w £ C°, £/ien a(w, it) = /3(iu, it) for every u £ R™ 
wii/i < F(x). 

Proof. Obviously, we have f3(w,u) a(w,u). For u e R" with F(u) < F(x), 
wc get 

a(w, u) = sup {w T y\ y £ -F(x) < -F(u)} 

< sup {w T y| y £ F(x), -F(x) < F(x)} = a(w, x). 

Since F(u) < F{x) can be written as F(u) 2 -F(x) + C and, as id e C°, we 
obtain 

/3(tu, u) = sup {iu T y| y £ -F'(u)} 5= sup {u> T y| y £ -FXx)} = /3(u>, x). 

The point x being a solution of ( |P(tt;,x)| implies that there exists y £ F(x) 
such that P(vu,x) ^ w T y = a(w,x) ^ a(w,x). Altogether wc get a(w, it) ^ 
a(w,x) sg /3{w,x) ^ /3{w,u) si a(ui, u), which yields the desired equality. □ 

Theorem 4.4. SetOpt computes a solution of (P) whenever fP]) is feasible 
and bounded. Otherwise SetOpt states whether (|P]) is infeasible or unbounded. 

If the H-representation computed in phase 2 contains no redundant inequal- 
ities, SetOpt terminates after finitely many steps. To be more precise, let 
the pre-solution computed in the first phase consist of k points and let I be 
the number of linear inequalities necessary to describe the polyhedral convex set 
gr F + (Or" x C) . In the second phase of SetOpt, at most I ■ k linear programs 
have to be solved. 

Proof. If jP]) is infeasible or unbounded, so is (|VR|) and the algorithm terminates 
with the corresponding status. If ([P]) is feasible and bounded, a pre-solution 

{x 1 , x h }, k 5? 1 of (|P|) is obtained by solving (|VR[) . 

For fixed i £ {1, . . . , k}, denote by x 1 the value of the variable x % before the 
while loop has been entered and let x l be the value after the while loop has 
been left. Then, X := {x 1 ,...,x fc } is the pre-solution computed in phase 1 
and X := {x 1 , . . . , x fc } is the result of the algorithm. We will show that for all 
is {l,...,fc}, 

F(x l ) < F(x l ) and $ x e R" : F(x) S F(x l ) . (9) 
The first condition in (|9|) implies 

inf F(x) < inf F(x). 

xeX xgX 
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Since X is a pre-solution of ([P]), the infimum is attained in X, that is, 

inf F(x) = inf F(x). 

It follows that the infimum is also attained in X. The second condition in ((9]) 
states that X consists of only minimizers, whence X is a solution to (|Pj) . 

To show ©, let i e be fixed. The first condition of © fol- 

lows directly from the constraint F(x) < F(x l ) of the equivalent formulation 
(P(w 3 , x 1 )) of the linear program (P(w J ; y 1 , . . .y s ))- Note that by Lemma FI721 
an optimal solution of (P(w 3 ; y 1 ,.. -y s )) always exists in the case where ([P]) 
is feasible and bounded. The while loop is left only after flag has been set to 
zero. This requires r iterations in the inner for loop, which occurs only if for 
every j e {1 . . . , r}, a = a(w J , x 1 ) equals j3 = (3(w 3 , x l ). In case of w J /\\w 3 \\ e K . 
this is known by Lemma |4~51 i.e., (P(w 3 ; y 1 , . . . y s )) does not need to be solved. 
But a(w 3 ,x l ) = (3(w 3 ,x l ) for all j e {1 . . . , r} implies that x l is a minimizcr, 
compare Lemma 14. II 

To show finiteness, note first that there is a finite algorithm (such as Benson's 
algorithm) to solve (|VR[) in phase 1. Consider the map F : M™ =J W , F(x) : = 
F(x) + C. Of course, grF = grF + (O^n x C) is a polyhedral convex set. Thus 
it can be expressed as 



grF = Ux,y) el"x R q \ Ax + By ^ bj 
for some A e R ixn , B e R Lxq , b G M'. For every x e domF, we have 
F{x) + C = \y eW\ By > b - Ix] 



Consequently, every H-represcntation of + C that contains no redundant 
inequalities consists of at most / inequalities. With other words, there are at 
most I different vectors w 3 / \\w 3 1| for j e {1, . . . , k} in the algorithm, which proves 
the claim. □ 

The theorem immediately implies the following existence result. 

Corollary 4.5. If fP]) is feasible and bounded, a solution exists. 

To reduce the computational effort of the algorithm for specific problems, 
we suggest an additional rule. Let {(xi, y±), . . . , (xk, yk)} denote the solution of 
(jVR.p obtained in phase 1 and consider iteration i e {1, . . . , k] of the outer for 
loop in phase 2: 

If yj £ F(xi) for j with i < j < k then skip all commands in iteration j. 

Clearly, this rule maintains the attainment of the infimum and thus the algo- 
rithm still works correctly. 

Finally we consider the special situation where F : ffi™ — » Qc, i.e., we have 
F(x) = F(x) + C for all ieP . In this case, © can be replaced by 

F(x) < F(x) F(x) 3 F(x) 

Vi = l,...,s : y l e F{x) 
Vi = l,...,s : Ax ^ b-By l 
<^^> > max \b — By l \ i = 1, . . . , s\ . 
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This means that the linear program (P(w; y 1 , . . . , y s ) I has only n + q variables 



and 2m constraints. The problem to obtain an H-representation of gr (F(-) + C) 
from an H-representation of gr F seems to be difficult in practice where it is 
typical that n » q, compare also Remark 12.11 One way to obtain it is vertex 
enumeration of a polyhedral convex set in K™xR'. In contrast, SetOpt involves 
vertex enumeration only in W . 
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